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High frequency dispersion does not alter the low frequency 
spectrum of Hawking radiation from a single black hole hori- 
zon, whether the dispersion entails subluminal or superlumi- 
nal group velocities. We show here that in the presence of 
an inner horizon as well as an outer horizon the superlumi- 
nal case differs dramatically however. The negative energy 
partners of Hawking quanta return to the outer horizon and 
stimulate more Hawking radiation if the field is bosonic or 
suppress it if the field is fermionic. This process leads to 
exponential growth or damping of the radiated flux and cor- 
relations among the quanta emitted at different times, unlike 
in the usual Hawking effect. These phenomena may be ob- 
servable in condensed matter black hole analogs that exhibit 
"superluminal" dispersion. 
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I. INTRODUCTION 

Recent work has shown that Hawking radiation is 
highly insensitive to modifications of the short distance 
physics of the quantum field. In these models linear fields 
are considered, and the field equation is modified at high 
wavevectors in some preferred frame, yielding a nonlinear 
dispersion relation u{k) relating frequency to wavevector. 
Models with both subluminal and superluminal [^|||] 
group velocities at high wavevectors have been studied, 
including lattice black hole spacetimes [|j (which have 
subluminal dispersion). The picture that emerges from 
these studies is that the thermal Hawking spectrum is 
very robust for black holes with temperature much less 
than the energy scale of the new physics. Although the 
short distance physics does modify this spectrum, the 
modifications are so slight at the frequencies of interest 
that they seem well nigh impossible to observe. 

We have found a dramatic exception to this rule how- 
ever. If there is both an outer and an inner horizon, and if 
the dispersion is superluminal, then the Hawking process 
for a bosonic field is self-amplifying and the radiated flux 
grows exponentially in time, while for a fermionic field 
the process is self-attenuating. What happens is that 
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the negative energy partner of a Hawking particle, af- 
ter falling to the inner horizon, "bounces" and returns 
to the outer horizon on a superluminal trajectory, where 
it either stimulates or suppresses more Hawking radia- 
tion in the bosonic or fermionic case respectively. This 
secondary radiation is not only different than the usual 
Hawking flux, but it is correlated to the prior radiation. 
In the bosonic case the process continues to amplify at 
least until the back reaction becomes important. 

Charged black holes have inner horizons, but astro- 
physical ones would loose their charge very rapidly, so it 
is difficult to imagine how this runaway Hawking effect 
could ever be observed for real black holes. Even so, it 
provides an interesting theoretical laboratory in which to 
explore the effects of short distance physics. Moreover, 
it is conceivably relevant to string theory, and it might 
be observable in a condensed matter analogue of a black 
hole. Let us briefly indicate these ideas in turn. 

In spite of many points of close agreement between 
the physics of near extremal D-branes and black holes, a 
glaring discrepancy persists. If a radiating near extremal 
D-brane state is maintained at fixed energy by a constant 
influx of energy in a pure state, then the entropy in the 
radiation will be constant and there will be correlations in 
the radiation that emerges at different times. For a black 
hole, on the other hand, the usual Hawking process leads 
to uncorrelated thermal radiation for all time. The effects 
of superluminal dispersion invalidate the usual Hawking 
picture because the negative energy partners return to 
the event horizon. If there is something analogous to the 
superluminal dispersion of our model in string theory, 
then perhaps that could eliminate the discrepancy be- 
tween the string and black hole pictures. This may not 
be so far-fetched. String theory is, after all, non-local 
in some sense, and there is some evidence Q suggesting 
that it supports superluminal effects. 

A condensed matter analog — Unruh's sonic black hole 
lll'0jl] — '^^s the original stimulus for the development of 
the dispersive models. In this model, a sonic horizon oc- 
curs where the flow velocity of an inhomogeneous fluid 
exceeds the speed of sound. Although it seems unlikely 
that this situation can be experimentally realized for a 
low temperature superfluid, there are variations of the 
idea that might be realizable, involving quasiparticles 
other than phonons in different systems. For example, 
this may occur for fermion quasiparticles in rotating su- 
perfluid vortex cores with gap nodes such as ■^He-A or 
d-wave superconductors [E0[, or in moving "^He-A tex- 
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tures In both these examples there are both inner 
and outer horizons. Moreover, the quasiparticle disper- 
sion relation is "relativistic" sufficiently near a gap node, 
and the group velocity increases (i.e. becomes "super- 
luminal" ) as the difference between the momentum and 
the gap node increases, so the effective field theory has 
"superluminal" dispersion. Thus it is not inconceivable 
that the phenomena discussed here may someday be ob- 
servable. 

This paper is organized as follows. In section || 
the superluminal dispersion model for both bosons and 
fermions is discussed. The propagation of wavepackcts in 
the black hole spacetime with in ner and outer horizons is 
analyzed qualitatively in section III and the implications 
for the amplification or suppression and the correlations 
in the Hawking radiation are drawn in section Section 
renders the previous discussion quantitative by using 
explicit wavepacket solutions (derived in the appendix) 
to find expressions for the number and correlations be- 
tween the radiated quanta. Open issues concerning the 
boundary conditions on the quantum state and the grav- 
itational back reaction are discussed briefly in section [v| . 

We use units with 7i — c = 1 and metric signature 
(+ )■ 



II. SUPERLUMINAL DISPERSION MODEL 

A 2-dimensional model suffices to illustrate the essen- 
tial physics. We assume the spacetime metric is static, 
and therefore jl2j coordinates can be chosen (at least lo- 
cally) so that the line element takes the form 



ds^ — dt^ — {dx — v{x)dt) 



(2.1) 



A special case is the line element of the t — r sub- 
space of the Reissner-Nordstrom black hole space- 
time in Painleve-GuUstrand coordinates, where v(r) = 
-v/2GM/r- QVr2. (These coordinates cover the black 
hole interior down to where v{r) = 0, at r = Q"^ /2GM .) 
More generally, we consider any v{x) which is negative, 
vanishes as x ^ +oo, and is greater than —1 except be- 
tween inner and outer horizons, located at Xi and Xo^ 
where v(xi^o) = — 1- 



A. Boson field 

We adopt a linear field theory with higher spatial 
derivatives included in the action in order to provide a 
superluminal dispersion relation. In this section we re- 
strict to the case of a real bosonic field. The case of a 



Majorana fermion field will be discussed in section II B 
The action for the field is given by 



Si, 



d^x 



(2.2) 



In the ordinary relativistic action one has F{dx) 
In this paper we take 

i?(a.) ^dl^ ^dt 



(2.3) 



To motivate this action we note that the black hole de- 
fines a preferred frame, the frame of freely falling ob- 
servers. In the Painleve-GuUstrand coordinate system, 
{dt + vdx ) is the unit tangent to free fall world lines that 
start from rest at infinity, and dx is its unit, outward 
pointing normal. Our action comes from modifying the 
derivative operator, only along the unit normal 0^, by the 
addition of higher derivative terms which become impor- 
tant only when the wavelength is of order 1/fco or shorter. 
We will assume that this length scale of "new physics" 
is much shorter than the length scale of the metric (^.l|), 
i.e. fco ^ k'/^l- (III particular, we assume fco S> n, 
where n — \v'{xi^o) \ is the surface gravity of the horizon.) 
The idea is that the microstructure of spacetime, or of a 
condensed matter analog, might give rise to such higher 
derivative terms in the effective action. The choice ( |2.3| ) 
is just the generic form for the lowest order such term that 
is reflection invariant and produces superluminal group 
velocities. 



The action (2.2,2.3) produces the equation of motion 



{dt + dxv){dt + vdx)4> = dl4> - -^dl<t>. (2.4) 

To derive the dispersion relation for this equation we look 
for solutions of the form 



(2.5) 



(/)(t, x) = exp —ibjt + i k{x') dx 



where k{x) is a position dependent wavevector. Substi- 
tuting this ansatz into the equation of motion (2.4) and 
neglecting derivatives of v{x) and k{x) results in the dis- 
persion relation 



where 



{uj-vkf ^F'^ik) 



F^{k)^e + k^/kl. 



(2.6) 



(2.7) 



The group velocity in the free-fall frame is dF/dk, so 
wavepackets with fc ^ fco propagate near the speed of 
light, whereas wavepackets with k }^ kg propagate super- 

luminally. 

The dispersion relation ( ^.6[ ) is a fourth order polyno- 
mial equation in the wavevector fc so it has four solu- 
tions for fc at given values of co and v. The nature of 
these roots is revealed by a graphical method. In flg- 
ure |l| we plot the straight line (w -I- |w|fc) for one value 
of Lo (satisfying < w <C fco) and two values of v, and 
the curve ±_F(fc), as functions of fc. (We define F{k) 
as the positive square root of (2.7).) The intersection 
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points are the allowed real wavevector roots to the dis- 
persion relation. When |d| < 1 there are only two real 
roots (corresponding to the two roots for the ordinary 
dispersion relation with Ford{k) — ±/c), the other two 
being complex. The positive wavevector is denoted /c+s. 
When |w| > 1 + |(a;/fco)^/^ w 1, on the other hand, aU 
four roots are real, with one positive and three nega- 
tive. The positive wavevector is denoted fc^ in this case 
and, in decreasing magnitude, the negative wavevectors 
are denoted k- and k-s respectively (the other negative 
wavevector corresponds to an ingoing wave that plays no 
role in this paper so we do not give it a name). These 
roots are labeled in figure |^. 
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FIG. 1. Plot of |«|fc) (for one value of u and two values 
of v) and F{k) as functions of k. The intersection points of 
the curv es ar e the allowed wavevector roots of the dispersion 
relation (2.6). 



The dispersion relation plot in figure || is also quite 
convenient for tracing the motion of wavepackets in the 
background spacetime. The coordinate group velocity 
Vg — dx/dt of a wavepacket centered on a given wavevec- 
tor is given by 



duj I I ± 

dk dk ' 



(2.8) 



where ±dF/dk is the group velocity in the free-fall frame. 
Thus at any wavevector Vg is just the slope of the ±F{k) 
curve minus the slope of the straight line [uj + \v\k). 
For all four types of wavevectors fc±s,± of interest to us, 
±dF/dk is positive, hence the sign of Vg is determined by 
which of the two slopes is larger, something that is easily 
read from the figure. For w > 0, the group velocity for 
k+s and k± is positive, whereas for k-s it is negative. 

When generalized to a complex scalar field, the action 
(2.2) is invariant under constant phase transformations 
of the field. This implies the existence of a conserved 



current j^. The integral of the time component j" over 
a spatial slice serves as a conserved inner product when 
evalu ated on complex solutions to the equation of motion 
( ^.4|) . For the metric (2T), this inner product takes the 
form 

U,g)^i f dx(^ndt+vd,)g-g{dt+vd,)r), (2.9) 



where f{t,x) and g{t,x) are solutions to (]2^ 

Tw o classes of complex solutions to the field equation 
( [2.4| ) are of special interest for quantization. The first are 
the positive free fall frequency wavepackets. They can be 
written as sums of solutions satisfying 



{dt + vd,)f{t,x)^-iuj'f{t,x) 



(2.10) 



where ui' > 0. The second are the positive Killing fre- 
quency wavepackets. These are sums of solutions of the 
form e~^'^^(p{x) where uj > 0. A positive free fall fre- 
quency wavepacket confined to a constant v{x) interval 
at one time necessarily has a positive norm under ( ^.S| ), 
as does a positive Killing frequency wavepacket confined 
to a region where v{x) — (where Killing frequency co- 
incides with free-fall frequency). Since the norm is con- 
served, it is positive at all times if it is at one time, even 
when the wavepacket does not remain in an interval of 
constant or vanishing v{x). Note that if the wavelength is 
small compared to the scale of variations of v(x), then a 
positive free-fall frequency wavepacket will have positive 
norm even if v(x) is not constant. 

To quantize the field we assume that (j){t, x) is a self- 
adjoint operator solution to the field equation that sat- 
isfies the canonical commutation relations. We define 
the annihilation operator a{f) associated to a normal- 
ized complex solution to the wave equation f{t,x) by 



a(/)^(/,0). 



(2.11) 



The commutation relations for the field operator are 
equivalent to the relations 



K/),at(.g)] = (/,5) 



(2.12) 



for all / and g. If f{t,x) is a positive norm solution, 
then a(/) behaves as an annihilation operator. If f(t, x) 
is a negative norm solution, f*{t,x) has positive norm, 
so a{f) = —«'!'(/*) behaves as a creation operator. 



B. Fermion field 



For simplicity we consider two-dimensional massless 
Majorana fermions. Following the conventions of Jl3[, 
the action in a general curved spacetime is given by|^ 



^In higher dimensions there would be a spin connection term 
as well. In two-dimensions it is easy to show that this term 
vanishes identically. 
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(2.13) 



(2.20) 



where = F'^e^^ and is the zweibein. We take the 
flat space gamma matrices as 



-i 

1 



i 

1 



Decomposing the spinor ip as 



(2.14) 



(2.15) 



and expanding the action in the metric (|]^) using the 
zweibein (ep, ei) = (pt + vdx, dx) we find 

^^ = \j d^^[-^+{dt + (1 + 

+V'-(at-(l-w)9,)V^_). (2.16) 

In this form it is clear that ^-i- and "0- do not mix. Fur- 
thermore at infinity, where v{x) = 0, V'-i- is right-moving 
while 'ijj- is left-moving. We therefore drop ip- in the 
remainder as it plays no role in the Hawking radiation 
calculation. 

Following the same motivation described in section ||, 
we now modify the action for by subtracting the 
higher derivative term k^'^ij^j^dxi^^^ obtaining the action 

S^^'-j d'x (V'+(5t + (1 + v)dx - K^dlyj+) . (2.17) 

Varying with respect to results in the equation of 
motion 

(dt + vdx + dxv/2 + dx- fco"'9^) V'+ = 0. (2.18) 

Substituting ip+{t,x) — exp{—iijjt + i k{x')dx') into 
the equation of motion and dropping derivatives of k{x) 
and v{x) results in the dispersion relation 



vk 



k^ /kQ. 



(2.19) 



This is the same (up to the coefficient of the k^ term 
and higher order terms) as the branch of the scalar field 
dispersion relation corresponding to pos itive group ve- 
locity in the free-fall frame given in ( p.q ) and displayed 
in figure |^. The classification of scalar wavepacket types 
in section ( |ll A| ) therefore applies to fermion wavepackets 
as well. In particular, the higher derivative term leads to 
super luminal propagation at large wavevectors. 

To quantize the field we assume that ip^{t,x) is 
a self-adjoint operator solution to the field equation 
that satisfies the canonical anti-commutation relations 
{'ilj+{t,x),'tp+{t,x')} = 6{x,x')}. The conserved inner 
product is the integral of the time component of the con- 
served current associated with phase invariance of the ac- 
tion (2.17) (generalized to complex fermions), and takes 
the form 



We define the annihilation operator b{f) associated to a 
normalized complex solution to the wave equation f{t, x) 
by 



(2.21) 



The anti-commutation relations for the field operator are 
then equivalent to the relations 



{b{f),b^{g)}^{f,g) 



(2.22) 



for all / and g. We represent the operators 6(/) on the 
fermionic Fock space generated by positive fr ee-fal l fre- 
quency solutions to the equation of motion (2.18). If 
f{t,x) is a positive free-fall frequency solution then b{f) 
behaves as an annihilation operator on this space. If 
f{t, x) is a negative free-fall frequency solution, then 
f*{t,x) has positive free-fall frequency, so b{f) = 6^(/*) 
behaves as a creation operator. 



III. WAVEPACKET PROPAGATION 

In this section we give a qualitative analysis of the 
role of the inner horizon in modifying the Hawking radi- 
ation. This analysis will exploit a WKB description of 
wavepacket propagation, allowing for non-WKB "mode 
conversion" in the vicinity of the horizons. The anal- 
ysis applies equally well for the bosonic and fermionic 
quantum fields. Scattering of waves on account of the 
background curvature of the metric ( p.l[ ) is negligible as 
long as the radius of curvature is much greater than 1 /kg. 
For small wavevectors, fc <|; fco, this is because the wave 
equation is approximately conformally invariant and the 
metric (hke any two-dimensional metric) is conformally 
fiat. For large wavevectors, k > fcoi it is because the 
wavelength is much smaller than the radius of curvature. 

We begin outside the outer horizon with a low fre- 
quency outgoing wavepacket peaked around a wavevector 
of type fc+s (see figure and we follow this wavepacket 
backwards in time. A sketch of what we find is given 
in figure ^. The final wavepacket (i.e. the one we begin 
with) is labeled +s in figure ^. This packet has posi- 
tive group velocity and therefore is right-moving, as can 
be seen from the graph of the dispersion relation (Fig. 

Following this packet backward in time it moves 
toward the black hole and blueshifts. The Killing fre- 
quency Lo is conserved, so the increase in the wavevec- 
tor can be seen from Fig. ^ by increasing the slope of 
the straight line while keeping the intercept fixed. As 
the wavevector grows, the group velocity increases in the 
free-faU frame, and so the packet becomes superluminal 
and crosses the horizon (backward in time), becoming a 
packet with wavevectors of type k+ (see figure 0). 
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FIG. 2. Spacetime sketch of the evolution of an outgoing 
k+s wavepacket backward in time. A line end indicates a 
wavepacket arising from mode conversion, while an unbroken 
line indicates that the wavevector evolves continuously on the 
dispersion curve. 

The wavepacket inside the horizon also has a fc_ com- 
ponent, which is not obvious if we simply follow contin- 
uously along the dispersion curve. In fact, the WKB ap- 
proximation breaks down near the horizon, and "mode 
conversion" from the positive wavevector to the nega- 
tive wavevector, negative free-fall frequency, branch of 
the dispersion relation occurs. This is easily shown ana- 
lytically, and is made plausible by the fact that, around 
the horizon, the straight line of figure ^ nearly coincides 
with a large portion of the curved line of the dispersion 
curve, thus allowing other wavevectors to become mixed 
in. The dispersion relation allows wavevectors of types 
fc_, and k-s in between the horizons, however only 
the first two are right moving, whereas the last type is 
left moving. Since our final wavepacket is by assumption 
purely outgoing outside the horizon, there can be no k-s 
component generated here. The /c+ and fc_ wavepackets 
are labeled -|- and — in figure I In this fi gure a line 
end indicates a wavepacket arising from mode conver- 
sion, while an unbroken line indicates that the wavevec- 
tor evolves continuously on the dispersion curve. 

The A:+ and fc_ packets propagate backward in time to- 
ward the inner horizon where they both undergo partial 
mode conversion. The group velocity of the fc+ packet re- 
mains positive around the inner horizon and therefore it 
can cross, becoming a fc+s packet, labeled +S2 in figure 
^. As before though, there is also some mode conver- 
sion from the positive to the negative wavevector branch 
of the dispersion relation, and a left-moving k^s packet 
(— S2 in figure 1^) is generated which propagates backward 
in time back toward the outer horizon. The fc_ packet 
on the other hand cannot cross the inner horizon on the 



negative wavevector branch because its group velocity 
drops to zero at the horizon. Indeed the group velocity 
goes through zero and becomes negative, so the /c_ packet 
turns around and propagates back toward the outer hori- 
zon as a k-s packet still on the negative wavevector 
branch. In addition, some mode conversion from the neg- 
ative to the positive wavevector branch of the dispersion 
relation occurs at the inner horizon. Therefore part of 
the fc_ packet does cross the horizon as a fc+s packet and 
is superposed with the fc+s packet that evolved from the 
k+ packet. 

The k+a packet inside the inner horizon continues 
propagating to the left backward in time. The k-s packet 
however returns to the outer horizon, near which its 
group velocity drops to zero. Again, partial mode con- 
version to the positive wavevector branch occurs, so the 
/c_s packet evolves backward in time to a pair of fc+ and 
/c_ packets which are heading back to the inner horizon. 
This is now almost the same situation we started with, 
since the original /c+s packet also evolved into a pair of 
k^ and fc_ packets between the horizons (although with a 
different relative weight). The analysis given above thus 
tells us qualitatively what happens when they reach the 
inner horizon, namely, the same thing as happened be- 
fore. The general pattern that emerges is shown in figure 

I 

We have so far discussed the history of an outgoing /c+s 
wavepacket followed backward in time. It is also instruc- 
tive to look at the future evolution of a k-g wavepacket 
in between the horizons, since the negative energy part- 
ner of a Hawking particle is such a wavepacket. This 
evolution can be inferred by the same sort of analysis 
just given, or simply by time and space reversal of that 
analysis, and is shown in figure |^. 
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FIG. 3. Spacetime sketch of the trajectory of a Hawking 
particle and its partner forward in time. A hne end indicates a 
wavepacket arising from mode conversion, while an unbroken 
line indicates that the wavevector evolves continuously on the 
dispersion curve. 



IV. PARTICLE CREATION: ORIGIN OF THE 
AMPLIFICATION OR SUPPRESSION OF 
HAWKING RADIATION 

The amount of particle creation in an outgoing posi- 
tive frequency wavepacket tp is indicated by the expec- 
tation value of the number operator N{'ijj) = {ip)a{il;) . 
To determine this expectation value an initial quantum 
state must be specified. Let us define an in-Hilbert space 
on some spacelike surface as the Fock space generated 
by positive free-fall frequency wavepackets on that sur- 
face. The corresponding ground state is then annihi- 
lated by annihilation operators of these wavepackets. We 
shall suppose the initial state is such a free-fall ground 
state associated with a given surface S. Decomposing 
ijj = ip^ + ip~ into its positive and negative free-fall fre- 
quency parts on E, the ground state condition implies in 
the bosonic case that {N{ip)) — and in the 

fermionic case (7V('0)) = {ip~,'ip~). 

Suppose we choose S as surface 1 in figure ^, i.e., a 
surface that cuts through the fc+ and fc_ packets first 
produced by propagating the fc+s packet back in time. 
Then the number expectation value for the fc+s packet 
is just (minus) the norm of the fc_ packet. In [^Jg] this 
was shown (for bosons) to be thermal at the Hawking 
temperature, for wavepackets with Killing frequencies m 
satisfying k ^ w ^ fcg. That is, the standard Hawking 
effect occurs even in the presence of superluminal disper- 
sion, if there is only a single horizon. 

When there is also an inner horizon, the particle cre- 
ation depends very much on which surface is used to de- 
fine the initial ground state. If we impose the ground 
state condition on the earlier surface 2 in figure ^, in- 
stead of surface 1, the occupation number for the final 
fc+s packet is no longer thermal. The norm of the neg- 
ative frequency part of the wavepacket on surface 2 is 
determined not just by the final passage across the outer 
horizon, but also by the mode conversion processes at the 
inner and outer horizons. 

As the time between the initial ground state and the 
final outgoing wavepacket grows, there is an exponential 
amplification or suppression in the occupation number 
of the final wavepacket in the boson and fermion cases 
respectively. To see why, note that the k-s packet de- 
noted — S2 in figure Revolves into the orthogonal fc+s and 
k^s packets denoted -t-si and 
norms are related by 



-si respectively, hence the 



-S2\ 



-Si 



-Si 



(4.1) 



Consider first the bosonic case. A fc+s(fc_s) packet 
has positive (negative) free fall frequency and therefore 
positive (negative) norm under (p^), so it follows from 



(4.1) that II —Sill is larger in magnitude than || — S2I 



Continuing into the past this process repeats, and for 
each "bounce" between the horizons the norm of the 
wavepacket between the horizons grows by some fixed 
multiple, resulting in exponential growth of both 1 1 — s„ | p 



and II -I- 



Since the negative frequency part of this 



wavepacket determines the number of created particles in 
the final outgoing wavepacket, that number will grow ex- 
ponentially in the time between the initial surface S and 
the emergence of the outgoing wavepacket ij). Viewed 
forward in time, the Hawking effect is a self-amplifying 
process since the negative energy partners of the Hawk- 
ing particles return to the event horizon (in the form of 
a pair of fc+ and packets) and stimulate the emission 
of more radiation and more partners. The wavepacket 
trajectories associated with this forward in time picture 
are shown in figure |^. 

For a fermionic field, the above discussion is modified 
only by the fact that all wavepackets have positive norm. 



so equation (4.1) implies that || — si|| is smaller in mag- 



nitude than II — S2II . This means that the number of 
created particles will be exponentially damped in time. 
In effect, the allowed states between the horizons for the 
negative energy partners of the Hawking particles become 
filled, cutting off further pair creation. 

One further important point can be extracted from this 
analysis. Since a single particle/partner wavepacket pair 
evolves to a sequence of outgoing wavepackets as shown 
in figure ||, the states of all these outgoing wavepackets 
will be correlated. This can also be seen from the back- 
wards in time picture. It is clear from figure ^ that suc- 
cessive outgoing wavepackets will have past histories that 
partly overlap, in particular on the initial ground state 
surface, so there will be correlations between the quanta 
emitted from the horizon at different times. These corre- 
lations are in sharp contrast to the usual Hawking effect 
which produces uncorrelated thermal radiation. The in- 
formation loss that is normally associated with the cor- 
relations between Hawking quanta and their partners is 
largely eliminated, since an unending sequence of Hawk- 
ing quanta is coherently correlated to the same partner 
degrees of freedom. 



■^It is perhaps surprising to have exponential growth in time 
when there are no imaginary frequency solutions to the disper- 



sion relation (2.6). There is no contradiction however, since 



where 1 1/| ^ stands for (/, /) in the bosonic case and (/, /) 
in the fermionic case. 



these exponentially growing wavepacket solutions cannot be 
Fourier transformed in time, so need not be expressible as 
superpositions of time-independent mode solutions. 
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V. QUANTITATIVE ANALYSIS 

The qualitative analysis of the previous section 
will now be sharpened by explicitly constructing the 
wavepacket solutions discussed there. This will allow us 
to quantify the amount of amplification, suppression, and 
correlation of the black hole radiation. In the first two 
subsections we treat only the bosonic case, and in the 
last subsection we discuss the fermionic case. 



A. Wavepacket solutions 

The basic idea applied here is to patch together lo- 
cal wavepacket solutions with the aid of "evolution for- 
mulae" . The derivation of these evolution formulae is 
discussed in the Appendix of this paper. They are de- 
rived using connection formulae for time-independent 
mode functions which are obtained by matching WKB 
solutions to near- horizon approximations. Forming 
wavepackets with the mode functions we then obtain the 
evolution formulae for the wavepackets. 

Evolution formulae are needed for two different bound- 
ary conditions at both the inner and outer horizons, cor- 
responding to the spacetime diagrams in figures 4a-d. 
Using the notation "/ g " to denote "/ evolves to g " 
(forward in time) , the evolution formulae about the outer 
horizon are: 



Xn,+ + Xn 



while about the inner horizon they are 

n,—s ^ Xn,-\- Xn,— 



(5.1a) 
(5.1b) 



(5.1c) 
(5.1d) 



where all packets have been left unnormalized in order 
to keep the formulae as simple as possible. The evolu- 
tion formulae given here are preferred for evolving pack- 
ets backwards in time. Following the same techniques 
described in the appendix evolution formulae more con- 
ducive to evolving wavepackets forward in time can be 
derived. The -I-,— ,-|-s and —s subscripts denote which 
type of wavevector the packet is peaked about. ■(/'«, +s lies 
outside the outer horizon, while ^i^n,+s and Xi,n,+s both 
lie inside the inner horizon. The subscript n is a sort 
of time variable. Translation of n by one unit has the 
effect of translating the wavepacket in time by a certain 
amount and also distorting the wavepacket. Note that 



in the evolution formula at the outer horizon (5.1b) n 
increases by one unit on the tpn+i.-s wavepacket. The 
evolution formulae ( p.la -d) are basically scattering so- 
lutions about a black hole event horizon with ( ^.la -b) 
corresponding to scattering the 4'n,+s and ipn,-s packets 
off the outer horizon backward in time (see figures and 
b respectively) and (5.1c-d) corresponding to scattering 



the resulting combinations of -I- and — packets off the in- 
ner horizon backward in time as well (see figures ^ and 
d respectively). 



a. 




b. 



n+l,-s 





i,n,+s 



n,-s 




n,-s 



FIG. 4. Spacetime sketches of the local wavepacket evolu- 
tions given by ( 5.1a -d) respectively. 



To construct the wavepacket solution with final data 
consisting of a fc+s packet outside the outer horizon we 
start with the local solution (5.1a) (figure ^). This 
clearly is not a global solution since ?An,± do not solve 
the equation of motion (2.4) about t he in ner horizon. 
The combination may be replaced by (5.1c) which does 
however. This results in the evolution formula 



(5.2) 



The evolution (5.2) i s al so not a global solution to 
the equation of motion (2.4) however since ipn. -s is not 
a solution about the outer horizon. Using ( Frbl) to 
evolve i}}n,-s about the outer horizon, followed by (5. Id) 
to evolve the resulting wavepackets ■i/;„_i.+ and tpn-i,- 
about the inner horizon, we obtain an evolution formula 
that can be iterated indefinitely, 



(Xi,n- 



1,-s 



) ^ V-n, 



(5.3) 



Beginning with (5.2) and iterating (5.3) (n — m) times 
yields 



Xi,n—j,+s + Wm.- 



+8- 



(5.4) 



In this manner we can evolve the final wavepacket packet 
back to the spacelike surface where the initial ground 
state boundary condition is defined. Solutions of the 
form (x4) correspond to those used in the qualitative 
discussion of section III, which are depicted in Fig. 
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(wherein t he firs t fc_s packet has been traded for a k± 
pair using ( 5.1b ).) 

If the horizons are not sufhciently widely separated 
then intermediate wavep acke ts that arise between the ini- 
tial and final packets of (5.4) will overlap with the initial 
and final wavepackets, thus complicating the analysis of 
particle creation. We can avoid such complications by 
constructing a different solution. Setting n equal to m in 
(5.2) and subtracting from (5.4), we obtain 



n—m 



Xi,n—j,+s Wi,ra,+s 

- V'm,+s)- (5.5) 



This solution corresponds to sending (n—m) Xi,k,+s pack- 
ets and a pair of ipi^k,+s packets into the inner horizon 
and getting a pair of V'+s packets out of the outer horizon. 



B. Particle creation 

We can now compute the average number of particles 
in the wavepacket^ '4'n,+s and the correlations between 
emitted particles for different values of the time index 
n. To begin with let us evaluate the occupation number 
(0|A^(V'o,-i-s)|0) of the first outgoing packet after the ini- 
tial state condition is assumed. Then the only evolution 
formula we need is (5.1a), with n = 0. The annihila- 
tion operator for a normalized wavepacket / is given by 
fl(/) =^ {!:'!>) ( 2.11 ). Taking the inner product of ( 5.1a ) 
with the quantum field 0, and using the ground state 
conditions 



a(Vo.+)|0) =0-a(V3oV)|0), 



(5.6) 



we obtain 



(0|iV(^o,+.)|0) = - 



(-00 -,-00 -) 
(■0ri,+s, 0n,+s) 

duj 



LUu - UJl 



g27raj/ K ^ ' 



(5.7) 



(The norms ( A30 ) were used in the last equality.) This is 
just the Planck distribution at the Hawking temperature 
Th = k/2tt, as was shown previously for a superu- 
minal dispersive field theory in the case that there is just 
one horizon. It holds for k < w <C fco- 

We would have obtained a different result for 
(Oj_/V(0o,+s)|O) had we replaced the initial conditions 
( p.6|) with, for example, 



a(^,,o,+.)|0) =0 = a(4*_J|0>• 



(5.8) 



^We use a 'to denote normalized wavepackets. 



Indeed, from (5.2) with n = 0, we find that if ( pT^ ) holds 
the occupation number of -0o,+s is given by 



(0|iV(^o,+.)|0) = 



(V^O,-^,0O,-s) 
(■00, + s, -00,-^5 ) 
1 



du! 



2(1 - cos(6'+(u;) - 0-iuj))) 



(5.9) 



The phase angle s 6*+ (ii;) are defined implicitly in ( Al^ ), 
and the norms ( [A30| ) were used in the la st e quality of 
( ^.91 ). This differs from the thermal result (^.7|). 

It is not yet clear to us what is the "correct" initial 
condition on the quantum state of the field. To deter- 
mine this would require following the evolution of the 
field state as the black hole (or condensed matter black 
hole analog) forms. It does seem however that the con- 
ditions a(-0i^fe^+i)|O) — are likely to hold, while the 
remaining specification of the state remains to be de- 
termined. Fortunately these conditions alone sufhce to 
determine the rate of growth of the number of particles 
emitted and the correlations between them. 

In order to find the number and correlations in the 
radiation for n > we use the solution (|]^). We take n > 
m ^ 1 so that the intermediate wavepackets that entered 
the construction of ( |5.5D will not have any support on the 
initial surface, and we assume the ground state conditions 



a{ipi^k,+s)\0) = = a{xi,k,+s)\0) 



(5.10) 



for k > m. Taking the inner product of (|5.5| ) with the 
field operator (j) and using conditions ( 5.10 ) we obtain 



a{ipn,+s)\0) = 



(0n,+s, 0n,-(-s) 

from which it follows that 



a{lpm, + s)\0) , (5.11) 



{0\a\i^k,+s)a{4'n,+s)\0) = 

^(■0fc,+s, -0fe,+s)('0„,+s, 0n,+s) 

In particular taking k = n we obtain 

(0m, +s, ■0m, +s) 



(0|7V(V>™,+,)|0). (5.12) 



(0|7V(V„,+,)|0) 



(■0„,+s,-0, 



(0|iV(V3m,+.)|0). 



n,+s ) 



(5.13) 



The norm of ^n,+s (A30) is given by 

, f 1 - cos(6l+ - 6'_ 

(1pn,+s,'(Pn.+s) = 47r / dlV [ 1 + ^ ^ — - 

Juii \ 2smh [ttuj/k) 



(5.14) 



where we have used ( A21). This decreases monotonically 
with n except for at most a discrete set of frequencies for 
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+ = 9-Juj) + 27rk for some integer k. ( For t hese 



which 0. 

frequencies Ti ( A21 ) vanishes, so according to ( |A20| ) the 
corresponding mode is a bound state trapped between 
the horizons.) Therefore the particle creation in V'n.+s 
increases monotonically with n, diverging as n ^ oo. In 
particular, if th e wav epackets are narrowly peaked about 
a frequency w, (5.13) yields 



= 1 



1 — cos( 



5-) 



2sinh (ttcj/k) 



(5.15) 



(0|jV(^„,+,)|0) 
(0|iV(^„,,+,)|0) 

which grows exponentially with n — m. 

A measure of the correlation between emitted particles 
is given by 



C{m, n) 



\a'' ilpm,+s)a{lpn,+s)\0) 



1/2 



1, 



,)|0)(0|7V(^„^+^)|0) 

(5.16) 

independent of the difference n — m. This should be 
contrasted with the correlation obtained when (j) satisfies 
the ordinary wave equation. 



C{m,n) = ('0„i,+s,'0„,+s), 



(5.17) 



which is nonvanishing only to the extent that these 
wavepackets are not orthogonal. As n—m grows the over- 
lap of these wavepackets and hence th e corr elation (5.17) 
goes to zero, whereas the correlation (5.16) remains. 

Finally let us estimate the time between the succes- 
sive particle emissions (see figure ||) , i. e, th e difference in 
times when successive ipn,+s packets (A22) cross a fixed 
X coordinate. The trajectory of the packets is given ap- 
proximately by the condition of stationary phase. 



t : 



^ Arg (^v4„ ( w ) (x , w ) 



(5.18) 



and therefore the time At between the n*^ and (n -|- 1)*'* 
packets crossing the coordinate x is given approximately 
by 



At 



duj 



Arg(r3 



e ^^/''e 



(5.19) 



where we have substituted for using (A23), and is 



given by ( [A2l| ). 

Using the results given in the appendix the lo- 
dependence of the phase factors 7(0;) and 6'±(w) can be 
computed. Rather than carrying out this calculation — 
which we can only do explicitly for any particular v(x) in 
some approximation anyway — let us make a rough esti- 
mate. The interval At is determined by the time it takes 
a wavepacket to "bounce" back and forth between the 
horizons. If v{x) -I- 1 is of order unity between the hori- 
zons, then using the group velocity of the k-s and k± 
waves one finds that this bounce time is of order a, the 
coordinate distance between the horizons. 



C. Fermionic case 

In this section we briefly describe the differences in the 
quantitative analysis of the fermion case. The derivation 
of the wavepacket solutions for fermions parallels that 
given for the scalar field in section VA and yields an 
evolution formula very similar to (5.5). The final ex- 
pression for the aver age v alue of the number operator is 
identical in form to ( 5.13 ), however the norm iV„^-|_s now 
increases monotonically with n so that the number ex- 
pectation value decreases exponentially in n. This is to 
be expected since unlike the scalar case, the conserved 
norm (/, /) ( 2.20| ) is positive definite. The effect of this 
crucial difference is that instead of exponential growth of 
Hawking radiation we now get exponential decay. 



VI. DISCUSSION 

We left the question of the "correct" initial condition 
unanswered. For a condensed matter black hole it should 
be straightforward to deduce this by following the state 
of the field as the "black hole" forms. It seems fairly 
clear that the V'j.+s wavepackets inside the inner horizon 
will be in their ground states. What is less clear is the 
state of the wavepackets "01 s-iid "01^ between the 
horizons. For a real black hole — if one wants to entertain 
the possibility of superluminal dispersion — the same may 
be true. The V'i,-i-s wavepackets inside the inner horizon 
arise from ingoing waves that scatter around or through 
the central singularity of a Reissner-Nordstrom black hole 
in the manner discussed in [ p^ . Since the region inside 
the inner horizon is static, it would seem plausible that 
these are in their ground state as well. 

Another issue we have not touched is that of the grav- 
itational back-reaction to the radiation studied here. In 
the bosonic case the exponential growth of the number of 
negative energy Hawking partners between the horizons 
would surely rapidly entail a strong gravitational reac- 
tion. In the fermionic case, the exponential suppression 
of radiation leads quickly to a state with no radiation 
at all. This is hard to reconcile with the usual picture 
in which Hawking radiation is a robust consequence of 
a general "well-behaved" state near the horizon. One 
would expect that although the negative energy states 
of the Hawking partners in the ergoregion between the 
horizons become filled, there is not all that much energy 
in these states (since the partners at late times are the 
same as the partners at earlier times due to the "bounce" 
between the horizons) so the back-reaction should be lim- 
ited. If so, then why doesn't the Hawking radiation con- 
tinue? The answer, it would seem, is that although the 
state is reasonably well-behaved in terms of energy den- 
sity, it has peculiar features in just those modes relevant 
to the Hawking effect. 
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APPENDIX: WAVEPACKET SOLUTIONS 

In t his appendix we explain how the evolution formulae 
( ^.la -d) for wavepacket solutions are derived with the 
help of the results of [^. We treat only the bosonic case, 
although the fermionic case is essentially identical. 

The wavepacket evolution formulae are inferred from 
conn ectio n formulae for mode solutions to the field equa- 
tion (2^) of the form 



u{t, x) 



(Al) 



where (j){x,uj) satisfies the ordinary differential equation 
(ODE) 

- (x) + (l - v'^{x))(t)"{x) + 2v{x){iuj - v' {x))(t)' [x) 
-iLu{iLu - v'{x))(l){x) = 0. (A2) 

In H] such solutions were constructed for a black hole 
spacetime with a single horizon. The basic technique 
used was to find approximate solutions to (|A2| ) using 
the WKB approximation away from the horizon, and to 
match these solutions across the horizon by comparing 
to the near horizon solution obtained by the method of 
Laplace transforms. 



1. Outer horizon connection formulae 

Assuming that the horizon is located at x = 0, and 
that the metric in the vicinity of the horizon is given by 



v{x) ~ —1 + KX, 



(A3) 



the analysis of |^ leads to the following two connection 
formulae: 

K{e'"^/^^(j)+ + e-""/2'^'0_ ) ^ (j)+s (A4a) 
+ A'(e-'^"/2K.^_^ _^ e'^^/^";/).) ^ 0, (A4b) 



where 



K — ((jj/2 sinh(7raj/K)) 



.•^^1/2 



(A5) 



The notation "0i(a;) 02(2;)" denotes that the approx- 
imate WKB solution (l>i{x) behind the horizon connects 



to the approximate WKB solution (j>2 (x) outside the hori- 
zon . T he modes 0±, (f>±s are approximate WKB solutions 
to ( [A^ ) and are given byQ 



^±ix)^C±{v{xf-l)-'/^ 

X exp i / ds k±{v{s),uj) 



(A6) 



(t)±s{x) C±sexp J ds fc±s(u(s), w) j , (A7) 
where the approximate WKB wavevectors are given by 



k± « ±ko\/v'^ - 1 4- ujv/{l - v^) 



(A8) 
(A9) 



provided we assume w <C fco and choose \x\,e ^ 
{uj /kf)Y/'^ / K. The WKB solutions are only valid for 

|a;| > /^-l/3fco'^''^ The coefficients C±,±s are neces- 
sary to match these WKB solutions to the near-horizon 
Laplace transform solutions. They can be determined b 
comparing the Laplace transform solutions given in 
to t he m atching formulae (A4a-b) with the WKB modes 
(A^,A7) evaluated in the small x limit We find that 
the coefficients are given by 



C± 



_ ,-(lTl)/2 



exp 



N|3/2 



X exp 



=Fi-\/2K/fco(fcoe)^ 
(-z|^ln(2«e)) (AlO) 



C±s — exp ( i— In(fcoe) + i i argf r(l + iu / n) 

\ K 4 V 



(All) 



2. Inner horizon connection formulae 



In the case of a black hole with both inner and outer 



horizons, the connection formulae (A4a-b) remain valid 



*We have changed notation slightly from that in [gj. We 
have added the lower limit of integration ±e to the integrals 
appearing in the exponents and consequently the coefficients 
C±,±s acquire some e dependence to compensate. Further- 
more a factor of i appearing in the matching formulae of [c] 
has been absorbed in </!>_ and the phase of A*' as defined in f ] 
has been absorbed into (jjia- We have also renamed the (j>-m 
solution in |^ as <^_s here. 

cannot be arbitrarily small however because the WKB 
approximation breaks down as a; — > 0. In it was shown 
that the WKB and Laplace transform approximate solutions 
are both valid when K'^^^k^^^^ <^ \x\ <g and therefore 
the matching can be done in this range. We also choose e to 
satisfy the same inequality. 
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locally about each horizon (after some slight modifica- 
tions to be discussed presently), but the solutions are no 
longer global. Assume that the outer horizon is located 
at X q = 0, with v(x) taking the same form as given by 
(A3), and that the inner horizon is at Xi = —a with v(x) 
near the inner horizon taking the form 



v(x) K, —1 — k(x + a). 



(A12) 



(We assume that the surface gravity of the inner hori- 
zon has the same magnitude as that of the outer horizon 
to simplify the results. There is no difficulty however 
in allowing the surface gravities to be different.) Then 
t he connection formulae (A4a-b) are valid for equation 
( [A^ ) locally about the outer horizon, where the notation 
"01 (x) ^ 02(2^)" now denotes that 0i(a;) is valid between 
the horizons and (/)2 {x) is valid outside the outer horizon. 
To find the "local" mode solutions about the inner 



horizon we reexpress the mode equation (A2) in terms 
of the new coordinate y := —(x + a). The resulti ng y - 
equation is the complex conjugate of the x-equation ( A2), 
with v(x) replaced by v(i/) := v{—y — a). We denote the 

WKB mode solutions to this y-equation by (l)± and (l)±s ■ 
Since we have chosen the surface gravities to have the 
same magnitude, v{y) takes the same form near th e in - 
ner horizon as v{x) does near the outer horizon ( [A3| ). 
Therefore the mode solutions near the inner horizon arc 
the complex conjugates of those about the outer horizon, 
with X replaced by y, so they satisfy the same connection 
formulae: 



^ — TTUJ /2k. 



+ K{e' 



-I- e 



(Al3a) 

^^I'^'^'^S) ^ 0. (A13b) 



The notation "(^i(y) (^2(2/)" now denotes that 0i(y) 
is valid between the horizons and 02 (y) is valid inside 
the inner horizon. We have included a subscript "i" in 
(/)i_+s(y) to make it clear that this solution is valid only 
inside the inner horizon^ 

The WKB solutions 0±.-s(?/(a;)) are in fact the same 
functions of x, up to w-dependent phases, as the WKB 
solutions 0±^_s(x) respectively. To see th is, note that 
the WKB modes (/)±__s(x) given in (A6,A7) all take the 
formg 

(p{x) =C{uj)f{v{x))exp(^i J dsk{v{s),uj)^ (A14) 

where C is x-independent and / and k are real functions 
of v{x). Since, as discussed above, the x and y equations 
are related by substituting v{x) — > v{y) — u(— y — a) 
and complex conjugating, the WKB modes (j)±-s{y) are 
given by 



(f>{y)=C*iLu)f{viy))e^pi^-iJ ds k{vis),Lo) j (A15) 

Using y = —x — a an d vjy ) — w(x), and changing the 
integration variable in (A15) to s' — —s — a, we obtain 



-X — a) = C*(aj)/(w(x)) exp [ i / ds' k{v{s'),uj) 

-a+e 



(A16) 



which differs from (A14) only by an w-dependent phase 
factor, i.e., 

0(-x-a)= { ^C{uj) J dsk{v{s),u)j^ (/)(x). 

(A17) 

We shall not attempt to compute these phase factors, but 
rather shall assume the generic form 



-(-X - a) 



,{-x - a) 



0±(x), 



(A18) 



where 9±, 7 also depend on u(x), ko and a (the coordinate 
distance between the horizons) but do not depend on e 
as long as it is chosen large enough so that the WKB 
approximation holds. 



Using the phase r elations ( AIS) the inner horizon 
connection formulae ( A13a , Al3b|) can be reexpressed in 
terms of the same linear com binations of (j)± appearing 
in the outer horizon formulae ( A4a A4b| ): 



(A19) 

X(e-""/2K^^ + e""/2'^0„) -f T30_, <^ Ti4>r^+s (A20) 
with 



-) 



T2 = uj-'^K^{e'"^/'^e-''^+ - e^'^'^/'^e"*^-) 



— ,,,-1 f('2„i7/ -7rw/K„-j( 



'+ __gW«e-*«-). (A21) 



3. Wavepacket evolution formulae 

We can now form wavepackets from the modes and 
use the mode connection formulae to obtain wavepacket 
evolution formuale. To keep the latter simple, the w- 
dependcnce of the coefficients in the connection formulae 
is built into the definition of the wavepackets as follows. 
Define 



^We have dropped the ±, — s subscripts on (f>{x), C, /, and k. 
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V'n.+s 
Xn-l,± 



cLu 



i?„ e 



^f^„_ie-*'^V.,+s, (A22) 



where the coefficients A . . .F depend on tj, and the mode 
functions depend on both uj and x. 

With these definitions the evolution formulae about 



the outer horizon ( 5.1a . 5.1b ) follow im mediately from 
the outer horizon connection formulae (A4a.A4b) pro- 
vided we choose C„± = exp(±7rw/2K) A„ and Dn,± = 
-K exp(=F7ru;/ 2K) B n+i- Similarly, the inner horizon evo- 



lution formula (5.1c) follows provided i3„ = e'^TiA„ and 



En = T2An, while (5. Id) requires Fn = —T\Bn+\ and 
Bn = —T^Bn+i- The solution (up to an undetermined 
overall constant) is given by 



An = i-nr" 

S„ = e'^ri(-T3)-" 

i^„ = T2(-r3)"" 

Fn = -e''^(ri)2(-T3)-"-i 



(A23) 



4. Norm of the wavepackets 



The wavepackets defined in ( A22| ) are not normalized. 
Their norms can be determined as follows. A generic one 
of these wavepackets has the form 



(A24) 



which has the norm (cf. (2.9)) 



4,l{dt+vd,W^-cj,'^{dt+vd., 



(A25) 



The norm is conserved under time evolution so it suffices 
to evaluate it at any one time. 

The key assumption we need in order to evaluate the 
norm is that at some time the wavepacket is confined to 



a constant v region. This is certainly the case for the 
+s wavepackets, since they are outgoing and eventually 
reach the asymptotic region. If the region between the 
horizons is large and has a large constant velocity region, 
then it may similarly hold for the —s and ± wavepack- 
ets as well. Alternatively, these wavepackets spend some 
time squeezed near the horizon, with wavelengths much 
smaller than the length scale over which v{x) changes 
(but, in the case of fc^, still much longer than k^"^ , so 
we can nevertheless use the small k approximation). If 
the wavepacket is contained in a constant v region then, 
for the purposes of evaluating the norm, we can imag- 
ine this region to extend to infi nity in both directions. 
The fixed oj mode equation ( |A2| ) in a constant v region 
has solutions <j)^ = Cu) exp{ikx), where uj ~ vk — ±F(k) 



wit h F{ k) given by (2.7). Matching to the WKB modes 
(|A^,|A7|) we see that |C^| {1, {v^ - l)-^/*} for the ±s 



modes and ± modes respectively. Thus we have 



dx(/.:,(/.^ = 2^{l,(z;2-l)-3/4} 



duj 



dk 



S{uj' — ui). 

(A26) 



Using (|A26| ) in (|A25D yields 



(V',^)=47r{l,(z;2_i)-3/2} / dcolG^l 



.±F 



duj 



dk 



(A27) 

Using the small and large k approximations for k±s and 
k± respectively, we find that 



{l,(^;2_l)-3/2} 



±F 



duj 



dk 



±{l,c.-i} 



respectively. Thus, finally. 



(V',V) =±4^ / duo\G^\'{l,oj-^}. 



(A28) 



(A29) 



With (|A29|) a nd the coefficients (|A23|) for the 



wav epackets ( A22) we obtain the norms needed in section 



(V'o,->o.-) = -47r / duj 



1 



(V'„,+s,i/'n,+s) = 47r / duj\T'i\ 



-2n 



(V'o.-s, V'o.-s) = -47r / duj\Ti\ 



-An I duj 

I ull 



2(l-cos(6'+(w)-6'_(t^))) 
(A30) 
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